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! Abstract 
(N . 

' The spectrum of gravitational waves that have been produced in infla- 

. tion is modified during cosmological transitions. Large drops in the number 

i of relativistic particles, like during the QCD transition or at e~^e~ annihila- 

^ ' tion, lead to steps in the spectrum of gravitational waves. We calculate the 

O ■ transfer function for the differential energy density of gravitational waves 

I i for a first-order and for a crossover QCD transition. 

• ^ : 1 Introduction 

><: 

■ Detecting a stochastic background of gravitational waves ^ would open a new 

window to the early Universe. Primordial gravitational waves are predicted to be 
generated during inflation P| and could be detected with upcoming cosmic mi- 
crowave background observations. Defects, like cosmic strings, produce stochastic 
gravitational waves as well P]. Inflation and defects predict an almost scale- 
invariant energy density per logarithmic frequency interval for the most interest- 
ing frequencies (~ 10~® Hz for pulsar timing, ~ 10~^ Hz for LISA, and ~ 100 Hz 
for LIGO) of gravitational waves. 

The aim of this paper is to study the evolution of primordial gravitational 
waves through transitions of the equation of state, especially the QCD transi- 
tion. A step in the gravitational wave spectrum comes from the large drop of 
the number of relativistic particles (by a factor ~ 3) during the QCD transi- 
tion. Since entropy is conserved, the growth rate of the Hubble radius is 
diminished during the transition. Thus, the rate at which modes cross into the 
horizon is changed during the transition and a step in the spectrum shows up at 
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frequencies of order of the Hubble rate at the transition. For the QCD transition 
this frequency is 10^^ Hz today and the step is a 30% correction. Other large 
drops in entropy density happen at e^e~ annihilation (which gives rise to a 20% 
correction in the energy spectrum) and at a GUT phase transition. The typical 
frequencies are 10"^'' Hz for e^e~ annihilation and 10^ Hz for the GUT transition. 

Similar steps in the differential spectrum have been studied for gravitational 
waves generated by cosmic strings These gravitational waves are generated 
on subhorizon scales when the cosmic strings decay. Their frequency is > a~^H, 
where a is the ratio between the typical size of a string loop and the Hubble 
radius at formation of the loop, a is smaller than 0.1 and might be as small 
as 10^^ @]. In this situation steps in the differential spectrum follow from the 
conservation of entropy for decoupled species (gravitons) during cosmological 
phase transitions 0. This interpretation applies to modes that have been inside 
the horizon long before the transition or that have been generated on subhorizon 
scales. However, for superhorizon modes entropy and energy of a gravitational 
wave are not defined. We therefore cannot rely on the conservation of entropy 
argument when dealing with gravitational waves from inflation. 

In Sec. 2 we briefly recapitulate the generation of primordial gravitational 
waves from inflation. The equation of state during the QCD transition is dis- 
cussed in Sec. 3. In Sec. 4 we explain the origin of the step in the spectrum 
and present results from numerical calculations, for a first order QCD transition 
and for a QCD crossover. A further discussion of the steps in the spectrum of 
primordial gravitational waves from the QCD transition and other cosmological 
transitions, especially e~^e~ annihilation, is presented in Sec. 5. 



2 Primordial spectrum of gravitational waves 

Let us recall the production of gravitational waves during inflation 0, |^. The line 
element of gravitational waves is given by ds^ = —dt"^ + a?{5ij + hij)dx''dx^ . hij 
is a transverse, traceless tensor, a denotes the scale factor. The spatial average 
{hij{x)h^^ {x + r)) = J jQ{kr)k^\hk\'^d\n k defines the power spectrum k is 

the comoving wavenumber. We define the rms amplitude h oi a. gravitational 
wave per logarithmic frequency interval: h = k^^'^\hk\- The linearized equation of 
motion for h{t) reads 

h + 3Hh + —h = 0, (1) 

where the differentiation is taken with respect to cosmic time t and H = a/a. 

During the quasi-de Sitter period gravitational waves are produced with al- 
most scale-invariant spectrum. For superhorizon scales, fcph = k/a <^ iJ, the 
slow roll approximation gives |0: 
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where e = —H/H'^ ^ 1 is the slow roll parameter and F{t) = 1 — (7e + ln2 — 
l)e + 0{e^)\ 7e + ln2 — 1 0.27. The amplitude of gravitational waves stays 
constant until the second horizon crossing and decays as 1/a thereafter, h ~ 
Cfc sin(A;?7 + 5fc)/a, where r] = J dt/a is conformal time. Ck and 6k are determined 
by matching to the superhorizon solution. 

For subhorizon modes, fcph ^ H, the energy density of gravitational waves can 
be defined. The space-time average of the energy-momentum tensor over several 
wavelengths gives Pg = —{Mp/32n){hijh^^). The energy density per logarithmic 
interval in k is related to the rms amplitude h: 



k'-^ - ^klJ^h^ . (3) 



dfc 327r"P''2' 



The factor 1/2 comes from the time average over several oscillations. The energy 
fraction in gravitational waves, per logarithmic interval in k, is defined by 

(4) 

where pc = SMpHq/ISh). Finally, the fractional energy density per logarithmic 
frequency interval reads: 

By ~ we indicate that we replaced the exact time dependence of h by ahc/cto- 
Modes that cross into the horizon [| in the matter dominated epoch have a^d O'O = 
(ifo/^ph)^ and modes that cross into the horizon in the radiation dominated 
regime have ahc/«o = Ho/[{l + ZcqY^'^kpY^]. From Eq. (^) we obtain 

2 ^ fHds\' \ {H,/k^,,f Ho < fcph < Ho{l + z.^^/^ 



n,ik)^ ^^Fie)^^^J {(l + ^Iq)'! Ho{l +\qy/' < k,^Z H,s • 

For comparison with experimental limits we use the frequency today, / = 
2vrA;ph(to)- In inflationary cosmology the strongest limit on Qg comes from the 
anisotropics of the cosmic microwave background This limit yields ^ghf^ < 
3 X IQ-^^Ho/fy in the frequency range Hq < f < SOHq [|, i. From Eq. (|) 
we have ^ghl^ < 7 x 10~^^/2;eq ~ 1 x 10~^^ for modes that crossed into the 
horizon during the radiation dominated epoch. Direct limits on ^gh^Q in the 
radiation era have been obtained by pulsar timing (< 4 x 10~^ for / ^ 4 x 10~^ 
to 4 X 10~^ Hz) ll^, and from big bang nucleosynthesis (/ figd(ln/) < 10~^ 
for / > -ffBBN ~ 10~^ Hz) [|T2|. Although, inflation predicts an extremely small 
amount of energy in gravitational waves at small scales, this is different for other 



^We define the moment of horizon crossing by the condition fcph — H 
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models of structure formation. Defects might generate a scale- invariant spectrum 
of gravitational waves with amplitudes that may be seen in pulsar timing residuals 
in the near future. String cosmology predicts a different energy spectrum, oc 
thus the BBN constraint, which does not constrain the inflationary scenarios 
at all, turns into a severe restriction. 



3 The cosmological QCD transition 

Lattice QCD results indicate that QCD makes a transition from a phase of free 
quarks and gluons to the phase of hadrons at a temperature T^^ ~ 150 MeV [Q. 
This implies a Hubble radius i?H ~ 10 km at the transition. Modes that crossed 
the horizon during the QCD transition have frequencies f\ ~ 10~^ Hz today. The 
order of the QCD transition is still a subject of debate. Lattice QCD with Wilson 
quarks indicates that the QCD transition is of first order for the physical values 
of the quark masses JlSl, whereas results with staggered quarks ||16| indicate a 
crossover for the physical quark masses. 

We consider both scenarios for the cosmological QCD transition, a first-order 
phase transition and a smooth crossover. A first-order QCD transition starts 
with a short period of supercooling and the nucleation of hadron bubbles. The 
supercooling is tiny, 6T/T ~ 10"'^, which implies that the entropy production is 
negligible, AS/S ~ 10~^ ||18|, |l^]. After enough hadron bubbles have been nucle- 



ated to reheat the Universe to the critical temperature, further bubble formation 
is suppressed and the bubbles grow adiabatically during the remaining 99% of 
the transition. For our purpose the cosmic fluid remains in thermal equilibrium 
throughout the phase transition. From recent lattice QCD results [0 for the 
surface tension and for latent heat we find that the typical bubble separation is 
~ IQ-^ Ru li. 

A simple model for a first-order transition is given by the bag model . In 
this model the entropy density jumps at the critical temperature 



s{T) 
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1 + —e{T - 

fi'a 



(7) 



where the effective number of relativistic helicity degrees of freedom before the 
transition are = 51.25 (2 quark flavors, photons, and lepton^) and gs. = 17.25 
(pions, photons, and leptons) after the transition; A^f = — ga,. At the high 
and low temperature phases can coexist. 

For a smooth crossover we assume for the entropy density ^1 



s{T) 



Is 



1 + 



lA^ 

2 g. 



1 + tanh 



AT 



^e^, n"^, and all i^s. 



(8) 
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We recover the bag model for AT — 0. Both models coincide at temperatures 
far away from T^. 

In both models the pressure is determined from p = J sdT and the energy 
density follows from the second law of thermodynamics, p = sT — p. The relation 
between the temperature and the scale factor follows from the conservation of 
entropy, 

dlna ~ ~ds/dT ' ^ ^ 

except for T = T^, in the case of a first-order phase transition. In the bag model 
T (X 1/a for T 7^ T,,. During a first-order QCD transition, i.e. T = T*, the 
pressure p(Tj,) = p^ is constant. The energy density p(a) is obtained from the 
first law of thermodynamics dp = —3{p+p^)da/ a. In the bag model the transition 
lasts a third of a Hubble time. For a detailed discussion of the cosmological QCD 



transition see Ref. 19 



4 Step in the differential energy spectrum 

Figure |I] shows the transfer function Qg{f)/Qg{f <^ /^) from a numerical inte- 
gration of Eq. (|ip . The typical frequency scale is 

/^^i.sef^-V 10-^Hz, (10) 

V17.25/ 150MeV ^ ' 

which corresponds to the mode that crosses the Hubble horizon at the end of the 
bag model QCD transition. Scales that cross into the horizon after the transition 
(l.h.s. of the figure) are unaffected, whereas modes that cross the horizon before 
the transition are damped by an additional factor ^0.7. The modification of the 
differential spectrum has been calculated for a first-order QCD transition (bag 
model), Eq. (|^, and a crossover QCD transition, Eq. (||), with AT/T^ = 0.3. In 
both cases the step extends over one decade in frequency. The detailed form of 
the step is almost independent from the order of the transition. To solve Eq. (|1|) 
numerically, it proofs useful to use the scale factor as evolution parameter for the 
first order transition and the temperature for the crossover transition, a and T 
are related by Eq. (H). 

The size of the step can be calculated analytically: Long before and after the 
transition, in the radiation dominated era, 

H oc , (11) 

from the Friedmann equation and the conservation of entropy. On the other 
hand, /cph oc 1/a, independently of the equation of state. Thus, the rate at which 
modes enter the Hubble horizon is changing during the transition. This change 
leads to the step in the transfer function. 
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Figure 1: The modification of tlie energy density, per logaritlimic frequency in- 
terval, for primordial gravitational waves from the QCD transition. 

For a fixed mode fca that crosses the horizon after the phase transition the 
amplitude h is constant during the transition, the transfer function is one. Now, 
consider a mode /cb that crosses the horizon before the transition. Its amplitude 
decays proportionally to flhc/a- From Eq. @ we find 

^g{kh) oc kl^al^ oc g^^'^ , (12) 

where we have used Eq. (|TT]) and k^i^ = H at horizon crossing. The constants of 
proportionality have been neglected in Eq. ([T2|), because they drop out from the 
transfer function below. Comparing the differential energy spectrum for modes 
that cross into the horizon before the transition and for modes that cross into 
the horizon after the transition gives the ratio 

^4^= f^)^. 0.696, (13) 

for the QCD transition, which coincides with the numerical integration in Fig. |l[ 
The size and position of the step in the logarithmic spectrum is independent of 
the order of the transition! 

The result (13) is in agreement with the entropy conservation of subhorizon 



gravitational waves. However, for superhorizon modes the entropy is not defined. 

Let us estimate the slope of the step in the transfer function. From the 
Friedmann equation and the 1st law of thermodynamics we obtain the Hubble 
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rate function of the scale factor {w = p/p) 



d\nH{a) = -^[l+w{a)]d\na . (14) 

For < w < 1/3 during the phase transition the Hubble expansion changes 
during the transition, i.e., H oc with P = (5{a) in the range (3/2,2). This 
holds true if the change in w is small, i.e. (dw/da)Aa <^1 + w. To estimate the 
exponent 

<S) ^ ^ (15) 

we use the relation din / = dln(i/a)hc = [—(1 + 3w)/2 dlnajhc and Eq. (||). We 
find 

l + 3Mahe) ' ^''^ 

which takes values between —2 and for w between and 1/3. From the numeri- 
cal integration (see Fig. |^) we find that the minimum of the exponent is n ~ —0.8 
for the bag model, which corresponds to an effective value Wes ~ 0.15. This is 
consistent with the values w{a) G (0.09, 1/3) that are taken during the first-order 
QCD transition. For the crossover transition the drop in w is smaller and thus 
the step is not as steep as in the first-order transition. 



5 Discussion 

In Fig. |l] we indicated the frequency range (~ 1 yr~^) in which limits on fig 



have been reported from pulsar timing residuals |Tl|]. The frequencies where the 
step of the QCD transition would be visible is of the order 0.3 month"^. For 
pulsar timing the power spectrum of gravitational waves is more relevant than 
the energy spectrum. The power spectrum is oc flg{f)f~^. Our results show that 
the power spectrum might deviate from the behavior over a whole decade in 
frequency. Depending on the effective value of w during the transition it might be 
as steep as /^^ (for w = 0). However, a realistic estimate for the QCD transition 
gives Wrain > 0.1. Thus, the spectral index of the power spectrum, n — 4, lies, 
according to Eq. (|16|), between —6.1 and —5. 

The 30% step in the differential spectrum from the QCD transition might be 
bigger than the tilt of the spectrum of gravitational waves from inflation (see 
for a computation of the tilt for various inflationary models), e.g. in chaotic 
inflation. 

As mentioned in the introduction, similar steps in the spectrum might occur 
at e+e~ annihilation, where /ann ~ 10~^° Hz is the typical frequency. Eq. (^) 
is modified for e~^e~ annihilation, because neutrinos and photons are decoupled. 
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Taking their difference in temperature into account \Ty/T^ = (4/11)^/^] we obtain 

^g(/ > /ann) 
^^g(/ < /ann) 

where = 2,gu = 5.25, and gt = 10.75. Again, this result is in agreement with 
the entropy conservation of subhorizon gravitational waves. 

For the electroweak transition the typical frequency (~ 6 x 10"'' Hz) lies in 
the frequency range of LISA, however the modification of the spectrum is tiny 
(the only particle that disappears from the radiation fluid is the Higgs particle). 

We have studied the evolution of gravitational waves through the QCD tran- 
sition in some detail. Although the detection of primordial gravitational waves 
is unlikely in the near future, it is exciting that the particle content in the early 
Universe is remembered in the spectrum of gravitational waves today. The spec- 
trum contains information about the relativistic degrees of freedom at a certain 
epoch of the early Universe and the transition temperatures of various cosmolog- 
ical transitions. On the other hand, if we assume to know the QCD transition 
temperature from lattice QCD, measuring the position of the step in the differ- 
ential energy spectrum would tell us whether the gravitational waves have been 
created during inflation or in a more recent epoch of the Universe ~ 10^^ Hz 
from inflation, > 10~^ Hz from cosmic strings). 
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